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ABSTRACT. We analyzed almost complete 18S rDNA sequences of 10 bovine Eimeria species, namely Eimeria alabamensis, E. auburnensis, 
E. bovis, E. hukidnonensis , E. canadensis, E. cylindrica, E. ellipsoidalis, E. suhspherica, E. wyomingensis and E. zuemii. Although these 
sequences showed intraspecific variation in 8 species, the sequences of each species were clustered in monophyletic groups in all species, 
except E. auburnensis. The sequences constituted 3 distinct clusters in a phylogenetic tree with relatively high bootstrap values; however, 
the members including each cluster shared no similarities in oocyst morphology. 
KEYWORDS: 18S tDNA, Cattle, Eimeria, oocysts, sequences. 

doi: 10.1292/jvms.l3-0273; J. Vet. Med. Sci. 75(11): 1427-1431,2013 



Bovine coccidiosis causes acute and chronic intestinal 
damage to the cattle and causes considerable economic 
losses to the cattle industry. More than 20 Eimeria species 
have been identified, and least 13 Eimeria species are con- 
firmed as causative agents of this disease [9]. Pathogenicity 
to the host differs among the different Eimeria spp., and the 
2 species, E. bovis and E. zuernii, are highly pathogenic and 
sometimes cause severe diarrhea in calves. Identification of 
Eimeria species has been largely based on the morphology 
of sporulated oocysts; however, only some species show 
morphological resemblance with each other. Therefore, 
molecular characterization is essential to accurately delimit 
species and infer phylogenic relationships among species. 
The small subunit ribosomal RNA (SSU rRNA) gene is a 
well-known marker for molecular characterization. Several 
studies have addressed the phylogenetic relationships based 
on 18S rDNA sequences and the morphology of Eimeria 
species, which parasitize domestic fowl [2], rodents [16] and 
rabbits [6]. However, the complete 18S rDNA sequences of 
only 2 bovine Eimeria species {E. bovis and E. alabamensis) 
have been analyzed and, to date, no study has reported the 
relationships between phylogeny of bovine Eimeria spp. on 
the basis of 18S rDNA sequences and similarities in their 
oocyst morphology. The purpose of this study was to analyze 
almost the complete 18S rDNA sequence of 10 bovine Ei- 
meria species and determine the phylogenetic relationships 
among these species. 

MATERIALS AND METHODS 

Eimeria oocysts: One hundred sixty eight fecal samples 
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were obtained from cattle in Iwate, Aichi, Hyogo and Na- 
gasaki prefectures in Japan. Samples were evaluated for 
Eimeria oocysts using the sucrose centrifiagal floatation 
technique. Eimeria oocysts were isolated from positive 
samples, sporulated in 2% potassium dichromate solution at 
room temperature and morphologically identified according 
to the procedure described by Levine and Ivens [10, 11]. 
Oocysts of 10 Eimeria species, namely E. alabamensis, E. 
auburnensis, E. bovis, E. bukidtiotiensis, E. canadensis, E. 
cylindrica, E. ellipsoidalis, E. subspherica, E. wyomingensis 
and is. zuernii, were detected, and a total of 119 oocysts were 
used in this study (Table 1). 

DNA extraction: Total DNA was extracted from each 
oocyst. Briefly, a single oocyst was placed into a PGR tube 
using a capillary micropipette and an inverted microscope, 
washed 5 times with resuspension in 100 jul of 0.05% Tween 
20 solution and centrifugation at 10,000 rpm for 1 min and 
treated with 100 /// of 0.5% sodium hypochlorite for 1 hr 
at 4°C. After washing, the oocyst was treated with 100 /jI 
saturated sodium chloride for 1 hr at 55°C [5]. After washing 
several times with 100 /il MilliQ water and centrifiigation at 
10,000 rpm for 1 min, total DNA was extracted by exposing 
the oocyst to 5 cycles of freezing for 3 min at -80°C and 
thawing for 1 min at 60°C. The volume of total DNA solu- 
tion was adjusted to 6.65 jul with MilliQ water. 

PCR, sequencing and phylogenic analysis: Two primer 
sets for nested PCR, TK2 and ets2 for the first PCR, and TKl 
and 18S-14 for the second PCR, were designed to amplify 
an approximately 1,800-bp fragment of 18S rDNA on the 
basis of the reference sequences of E. bovis (EBU77084), 
E. alabamensis (AF291427), E. grids (AB544361) and E. 
reichenowi (AB544355) (Table 2). The first PCR reaction 
was performed in a 25-/// reaction volume containing 6.65 
Hl of DNA template, 0.2 mM of each dNTP, 0.1 mM of each 
primer, 0.25 U of KOD FX neo DNA polymerase (Toyobo, 
Osaka, Japan) and the manufacturer-supplied reaction buf- 
fer. The thermal program was 94°C for 10 min, 20 cycles of 
98°C for 10 sec, 67°C for 3 min and 72°C for 10 min. The 
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Table 1 . Profile of Eimeria oocysts used in this study 



Species 


Oocyst 




Location 


Isolate code 


Numbers 


h. alaoamensis 


hai-isolatel '-I—IU ' lU 


10 


Hyogo 


E. auhrnensis 


Eau-isolatel-l~7 


7 


Iwate 




Eau-isolate2-l~2 


2 


Aichi 




Eau-isolate3-l'-'4 


4 


Nagasaki 




Eau-isolate4- 1 —2 


2 


Hyogo 


E. hovis 


Ebo-isolatel-l~4 


4 


Iwate 




Ebo-isolate2-l~2 


2 


Aichi 




P nr»-i*;nl ^ IP's - 1 / 


2 


Aichi 




Ebo-isolate4- 1 — 3 


3 


Aichi 




Ebo-isolate5- 1 


1 


Aichi 




Ebo-isolate6-l~4 


4 


Nagasaki 




Ebo-isolate7-l~2 


2 


Hyogo 


E. bitkidnofisnsis 


Ebu-isolate 1-1 — 4 


4 


Iwate 




Ebu-isolate2- 1 —3 


3 


Nagasaki 




Ebu-isolate3-l~2 


2 


Hyogo 




Ebu-isolate4-l~3 


3 


Hyogo 


E. canadendsis 


Eca-isolatel-l~10 


10 


Iwate 


E. cylindrica 


Ecy-isolatel-l~4 


4 


Iwate 




Ecy-isolate2-l~7 


7 


Nagasaki 




Ecy-isolate3-l~2 


2 


Nagasaki 


E. ellipsoidciris 


Ee-isolate 1-1 ~2 


2 


Iwate 




Ee-isolate2-l~2 


2 


Iwate 




Ee-isolate3-l 


1 


Nagasaki 




Ee-isolate4-l~5 


5 


Hyogo 


E. suhspherica 


Es-isolatel-l~3 


3 


Nagasaki 




Es-isolate2-l~7 


7 


Nagasaki 


E. wyomingensis 


Ew-isolatel-l~10 


10 


Hyogo 


E. zuernii 


Ez-isolatel-l~4 


4 


Iwate 




Ez-isolate2-l~2 


2 


Aichi 




Ez-isolate3-l~4 


4 


Aichi 




Ez-isolate4- 1 


1 


Nagasaki 


Total no. 




119 





a) "isolate" means a fecal sample, b) "l~nos." mean oocyst no. 



Table 2. 


Primers used in this study 


Primers 


Sequences 


TK2 


5'-GGT TGATCC TGC GAG TAG TC-3' 


ets2 


5'-AAT CCC AAT GAA CGC GAC TCA-3' 


TKl 


5'-AGTAGT CAT ATG CTT GTC TC-3' 


18S-14 


5'-ACG GAA ACC GTG TTA CGA CT-3' 


18S-9 


5'-ACA ATT GGA GGG CAA GTC T-3' 


18S-13 


5'-ATG CCT ACC TTG GCT TCT CC-3' 


Ebl-R 


5'-AGC CTG CTT GAA ACA CTC T-3' 


18S-12 


5'-GAA CCG TAG TTC CCG ATC CT-3' 



second PCR reaction was performed using the same reaction 
mix as the first PCR, except 1 /i/ of the first PCR reactions 
was used as template DNA with Go Taq DNA polymerase 
(Promega, Madison, WI, U.S.A.). The thermal program 
was 94°C for 2 min, 30 cycles of 98°C for 10 sec, 56°C 



for 10 sec and 72°C for 90 sec, and 72°C for 7 min. PCR 
amplicons were purified using the High Pure PCR Cleanup 
Micro Kit (Roche Diagnostic, Mannheim, Germany) and 
directly sequenced in both directions using the BigDye Ter- 
minator v3.1 Cycle Sequencing Kit (Applied Biosystems, 
Foster City, CA, U.S.A.) using the additional forward prim- 
ers 18S-9 and 18S-13 and the reverse primers Ebl-R and 
18S-12 (Table 2). The resultant sequencing ladders were 
read using the ABI PRISM 3500-Avant Genetic Analyzer 
(Applied Biosystems). 18S rDNA sequences were aligned 
using Clustal X version 2.0 software [7]. A neighbor-joining 
(NJ) phylogram was constructed using PAUP 4.0bl0 [14]. 
For phylogenetic tree construction, the sequence of Eimeria 
tenella(a.ccession no. AF026388) was used as an outgroup. 
Bootstrap analyses were performed using 1,000 replicates. 
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Table 3. Pairwise similarity (%) of 18S rDNA sequences among 10 bovine Eimeria species 



Species (bp) 1 2 3 4 5 6 7 89 10 



1 . E. alabamensis 


1765- 


-1769 


98.8- 


-100 
































2. E. aubmensis 


1764 




97.0- 


-97.4 


99.8- 


-100 




























3. E. bovis 


1763- 


-1764 


96.7- 


-97.2 


98.8- 


-99.1 


99.8- 


-100 
























4. E. bukidnonensis 


1768- 


-1769 


98.4- 


-99.0 


97.3- 


-97.6 


96.9- 


-97.3 


99.5- 


-100 




















5. E. canadendsis 


1763- 


-1764 


96.8- 


-96.9 


98.9- 


-99.0 


98.2- 


-98.5 


96.8- 


-97.0 


99.9- 


-100 
















6. E. cylindrica 


1763- 


-1764 


97.0- 


-97.4 


99.5- 


-99.9 


98.7- 


-99.1 


97.3- 


-97.5 


98.8- 


-98.9 


99.8- 


-100 












1. E. ellipsoidaris 


1761- 


-1763 


96.7- 


-97.0 


98.8- 


-98.9 


99.3- 


-99.7 


96.9- 


-97.2 


98.3- 


-98.4 


98.8- 


-98.9 


99.8- 


-100 








8. E. subspherica 


1773- 


-1774 


98.4- 


-98.9 


97.3- 


-97.6 


96.8- 


-97.2 


98.3- 


-98.6 


97.0- 


-97.3 


97.1- 


-97.3 


96.9- 


-97.1 


99.8- 


-100 




9. E. wyomingensis 


1764 




96.8- 


-97.2 


99.5- 


-99.7 


98.8- 


-99.0 


97.0- 


-97.2 


98.8 




99.5- 


-99.7 


98.7- 


-98.8 


97.1- 


-97.2 


100 


10. E. zuernii 


1763 




96.9- 


-97.2 


98.9- 


-99.0 


99.4- 


-99.7 


97.1- 


-97.3 


98.4 




98.8- 


-99.0 


99.4- 


-99.5 


97.1- 


-97.2 


98.9 100 



RESULTS 

The 18S rDNA sequences of 10 bovine Eimeria species 
were determined and deposited into the DNA Data Bank of 
Japan (DDBJ) with accession nos. AB769547-AB769665. 
The sequences ranged from 1,761 bp fori", ellipsoidalis to 
1,774 bp for E. subspherica and showed intraspecific varia- 
tion in 8 species, except for E. wyomingensis and E. zuernii 
(Table 3). The variants included 6 for 10 oocysts of is. alaba- 
mensis, 3 for 15 oocysts ofE. auburnensis, 2 for 18 oocysts 
of E. bovis, 4 for 12 oocysts of E. bukidnonensis, 3 for 10 
oocysts of is. canadensis, 5 for 13 oocysts of is. cylindrical, 
3 for 10 oocysts of is. ellipsoidalis and 6 for 10 oocysts of is. 
subspherica. Sequence variation was relatively low (<0.2%) 
in 6 species, except for is. alabamensis (<1.2%) and is. 
bukidnonensis (<0.5%). The oocysts of is. wyomingensis and 
is. zuernii, which showed no intraspecific variation in their 
18S sequences, were obtained from 1 and 4 fecal samples 
(isolates), respectively; moreover, the feces from which is. 
zuernii oocysts were isolated were derived from 3 distinct 
locations (Iwate, Aichi and Nagasaki) (Table 1). Oocysts of 
is. alabamensis and is. canadensis, which showed intraspe- 
cific variation in their sequences, were each obtained from 
single fecal samples. The 18S rDNA sequences showed pair- 
wise similarities in the range of 99.5-99.9% (is. auburnensis 
vs. is. cylindrica) and 96.7-97.0% (is. alabamensis vs. is. 
ellipsoidalis) (Table 3). 

In the phylogenetic tree, the sequences of 9 Eimeria spe- 
cies, except for that of is. auburnensis, were clustered in 
monophyletic groups (clades) formed for each species, and 
the clades of is. alabamensis and is. bukidnonensis consisted 
of phylogenetically distant variants (Fig. 1). Ten Eimeria 
species constituted 3 distinct clusters with relatively high 
bootstrap values: cluster A contained of is. auburnensis, E. 
canadensis, E. cylindrica and is. wyomingensis, cluster B 
contained of is. bovis, E. ellipsoidalis and is. zuernii, and 
cluster C contained of is. alabamensis, E. bukidnonensis and 
is. subspherica. 

DISCUSSION 

Few studies have examined the relationships between the 
phylogeny of Eimeria species based on 18S rDNA sequences 



and their biological traits, such as oocyst morphology. Phy- 
logenetic relationships among Eimeria spp., which parasit- 
ize domestic fowl, correlate with similarities based on their 
oocyst size and shape [2]. Rabbit Eimeria species formed 
a phylogenic tree that was closely linked to the presence/ 
absence of oocyst residuum and showed no relevant simi- 
larities to other morphological characteristics, such as oocyst 
size and shape [6]. Although more than 13 Eimeria species 
parasitizing cattle possess biological characteristics unique 
to each species, few studies have addressed the relationships 
between their biological traits and molecular phylogenies. In 
this study, we clarified the 18S rDNA sequences of 10 Eime- 
ria spp. detected from cattle in Japan and found that these 
species were clearly divided into 3 monophyletic groups, 
although the members clustered in each group shared no 
similarities in oocyst morphology. It is reported that 6 bo- 
vine Eimeria species are clustered in 3 distinct groups in 
a phylogenetic analysis based on nuclear ribosomal ITSl 
sequences, although is. ellipsoidalis belonged to a group dif- 
ferent from that clustered in this study [4]. Interestingly, one 
of the groups (cluster B in this study) included the highly 
pathogenic is. bovis and is. zuernii. Similarly, these 2 species 
also constituted the same monophyletic group in analysis 
based on ITSl sequences [4]. Furthermore, the 18S rDNA 
sequences indicated that the highly pathogenic is. tenella 
and is. necatrix in domestic fowl also belonged to identi- 
cal phylogenetic groups [2]. The 2 chicken Eimeria species 
shared common characteristics, and their second-generation 
meronts develop within the lamina propria of the lower 
intestinal tracts [9]. The biological traits common to both 
is. bovis and is. zuernii are to parasitize the regions of the 
lower intestinal tracts, such as the cecum and colon and to 
produce first-generation meronts of large size [9]; these traits 
suggest that the 2 species may have evolved from a common 
ancestor, although is. ellipsoidalis which belongs to the same 
clade does not possess these traits. 

The results of our study revealed intraspecific sequence 
diversity in most bovine Eimeria species, particularly cru- 
cial differences (98.8% and 99.5%) in is. alabamensis and 
is. bukidnonensis. Nuclear ribosomal DNA in eukaryotic 
genomes consists of multiple tandemly arranged repeats, and 
the copies are thought to undergo concerted evolution, be- 
cause of unequal exchange and gene conversion within and 
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62 



100 



59; 



57 



Eau-isolate 

2- 1, 2-2, 3-2, 

3- 3, 3-4 

■ Eau-isolate 3-1 

Eau-isolate 
1-1, 1-2, 1-3, 
1-4, 1-5, 1-6, 
1-7, 4-1,4-2 



57 



50 



50 



Ecy-isolate 
2-2, 2-3, 2-4, 
2-5,2-6, 3-1,3-2 

Ecy-isolate 2-1 



U Ecy-isolate 1-1, 1-3, 1-4 

63'— Ecy-isolate 2-7 

67 I Ecy-isolate 1-2 

Ew-isolate 
1-1, 1-2, 1-3, 
1-4, 1-5, 1-6, 
1-7, 1-8, 1-9, 
1-10 



100 



Eca-isolate 
1-1, 1-2, 1-3, 
1-4, 1-5, 1-6, 
1-7, 1-10 



87 



99 



59 



— Eca-isolate 1-8 

-Eca-isolate 1-9 

Ebo-isolate 1-3 

Ebo-isolate 

1- 1, 1-2, 1-4, 

2- 1,2-2, 3-1, 

3- 2, 4-1,4-2, 

4- 3, 5-1,6-1, 

6- 2, 6-3, 6-4, 

7- 1, 7-2 



Ee-isolate 3-1 

Ee-isolate 
4-1 , 4-2, 4-3, 
4-4, 4-5 
Ee-isolate 

1- 1, 1-2, 

2- 1 , 2-2 
Ez-isolate 

1-1, 1-2, 1-3, 
1-4, 2-1, 2-2, 
3-1 , 3-2, 3-3, 
3-4, 4-1, 



99 

61 ' 



98 



90 r 



90 



.501 



67 



.81 



-Eal-isolate 1-10 



Eal-isolate 
1-1, 1-2, 1-3, 



1-4, 1-9 



Eal-isolate 1-7 



56 



99 



Eal-isolate 1-6 

Eal-isolate 1-8 
Eal-isolate 1-5 



9£ 



99 



95 Ebu-isolate 

\ 3-1,3-2,4-2,4-3, 



Ebu-isolate 

1- 1, 1-2, 1-4, 

2- 1, 2-3, 4-1, 
Ebu-isolate 1-3 

_ . -T- Etu-isolate 2-2 
Es-isolate 1^1 

Es-isolate 1-3 
Es-isolate 1-2 
Es-isolate 2-2 
Es-isolate 
2-3, 2-4, 2-5, 2-6, 2-7 
— Es-isolate 2-1 



99 




71 



521 



) 



0.005 



E. auburnensis 



E. cylindrica 



E. wyomingensis 



E. canadensis 



E. bovis 



E. ellipsoidalis 

E. zuernii 

E. alabamensis 

E. bukidnonensis 
E. subspherica 

£. tene//a(AF026388) 



A 



B 



C 



Fig. 1. Neighbor-joining phylogenic tree based on 18S rDNA sequences of 10 bovine Eimeria species. Each isolate shows a sequence 
obtained from total DNA extracted from a single oocyst profiled in Table 1 . The numbers in the branches represent bootstrap values of 
more than 50%. A, B and C represent the monophyletic groups, including the sequences of Eimeria spp. 
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between the same chromosome loci and appear to be contin- 
uously homogenized [1, 3]. However, several apicomplexan 
species including Eimeria species shared distinct types of 
18S rDNA with sequence identities ranging from 89% for 
Plasmodium falciparum [12] to 99% for Cryptosporidium 
parvum [8], Babesia bigemina [13] wA Eimeria mitis [15]. 
Although the distinct types are thought to be intragenomic 
variants, we were unable to clarify whether intraspecific 
sequence diversity coexists in a single genome since total 
DNA extracted from a single oocyst produced unique 18S 
rDNA sequences. Further molecular studies are required to 
examine this possibility. 
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